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• We study field theoretical models for cosmic (p,q)-superstrings in a curved space-time. We discuss 

, both string solutions, i.e. solutions with a conical deficit, but also so-called Melvin solutions, which 

I have a completely different asymptotic behaviour. We show that globally regular gravitating (p,q)- 

^T*' strings exist only in a finite domain of the parameter space and study the dependence of the domain 

. of existence on the parameters in the model. We find that due to the interaction between strings, 

the parameter range where string solution exist is wider than for non-interacting strings. 

o ■ 

m ; 

^ • I. INTRODUCTION 

^ i The recent realization that some extended objects in fundamental string theory [1] can have cosmolog- 

I— I , ical size has triggered a lot of interest both in the fundamental string and the cosmic string community. 

I , In braneworld scenarios, big quantities of two types of cosmic size "strings" (cosmic superstrings) can be 

\ produced Q when bane-antibrane pairs annihilate each other: fundamental (F) strings and Dirichlet (D) 

' strings, which can be either a one dimensional brane (Dl-brane) or a higher dimensional brane with all 

' but one dimension compactified. These two kinds of strings interact with each other Q forming bound 

' states: p F-strings can interact with q D-strings to form (p,q) bound states, resulting in a network of F-, 

■ D- and FD-strings. 

00 , Networks of (p,q) string have been investigated in several papers lately, both analytically and numer- 

. . , ically, and both in the realms of fundamental strings and solitonic gauge strings BS, 9,0, B- The 

. ^ I dynamics of (p,q) network is more complicated than "usual" (solitonic) cosmic strings: The intercom- 

' mutation of solitonic strings is almost always one, whereas in (p,q) strings it is usually less than one, 

' and often much less than one. Generally, solitonic strings do not form bound states, so that 3-way Y 

junctions do not form, whereas (p,q) strings are known to form stable bound states, and therefore create 
Y junctions. All strings in a network of ordinary solitonic strings have the same tension, whereas there 
is a whole range of different tensions in a (p,q) network. 

Properties of ordinary (solitonic) comics strings have been thoroughly studied, and their network 
properties are better understood than those of (p,q) strings. They are known to b e g enerically formed 
in supersymmetric hybrid inflation [9] and grand unified based inflationary models [lO|. There has been 
considerable effort in numerically modelling networks to obtain observables; for example, modeling fleld 
theoretical cosmic strings to obtain CMB power spectra 
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The gravitational properties of solitonic cosmic strings have also be studied in detail by minimally 
coupling the Abelian Higgs model to gravity. Far away from the core of the string, the space-time has a 



deficit angle, i.e. corresponds to Minkowski space-time minus a wedge [1^. The deficit angle is in linear 
order proportional to the energy per unit length of the string. If the vacuum expectation value (vev) 
of the Higgs field is sufficiently large (corresponding to very heavy strings that have formed at energies 
much bigger than the GUT scale), the deficit angle becomes larger than 2tt. These solutions are the 
so-called "supermassive strings" studied in [13] and possess a singularity at a maximal value of the radial 
coordinate, at which the angular part of the metric vanishes. Interestingly, it was realized only a few 
years ago Q that both the globally regular string solution as well as the supermassive solution have 
"shadow" solutions that exist for the same parameter values. The string-like solutions with deficit angle 
< 2tt have a shadow solution in the form of so-called Melvin solutions, which have a different asymptotic 
behaviour of the metric and have higher energies than their string counterparts (and are thus very likely 
cosmologically not relevant). The supermassive string solutions on the other hand have shadow solutions 
of Kasner-type. Kasner solutions possess also a singularity at some finite radial distance, the difference 
is that for Kasner solutions, the ti-component of the metric vanishes at this maximal radial distance, 
while the angular part of the metric diverges there. Since both supermassive as well as Kasner solutions 
contain space-time singularities they are surely of limited interest for cosmological applications. 

It is partly due to the fact that advanced machinery exists to study solitonic strings that different 
authors have tried to capture properties of (p,q) string networks using field theoretical models. One can 
construct field theory models that give rise to solitonic strings that share some of the properties of cosmic 
superstrings. For example, Abelian Higgs strings can have a reconnexion probability which is less than 



one when the strings meet at very high velocities [16[. Also, type-I Abelian Higgs strings can form bound 
states (and have strings of different tensions) when they meet at low velocities 0, or if a network 
of (type I) strings with high winding numbers forms in a thermal phase transition 18]. In models with 
SUSY flat directions strings with type I properties can also be formed 0. 

A maybe more realistic approximation is given by the models of Saffin Q and Rajantie, Sakellari- 
adou and Stoica 0], where two independent Abelian Higgs string models are coupled via an interacting 
potential. The strength with which the strings interact and form bound states can be controlled by 
the parameters in the potential, giving rise to numerically more feasible and controllable simulations, 
including network simulations [a, |8| . 

In this paper we couple minimally the models of Saffin and Rajantie et al. to gravity, and investigate 
the gravitating properties of this field theoretical approximations of (p,q) cosmic superstrings. In the 
next section we describe the models and set up the ansatze and boundary conditions. In section Hill we 
report the numerical solutions found for the models studied, and then conclude. 

II. THE MODEL 

The models we are studying are given by the following action: 
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where R is the Ricci scalar and G denotes Newton's constant. The matter Lagrangian £,„ reads: 

£,„ - D,<PiD^^r - If^uF^'' + D^aD^'O* - -^H^.H^"" - V{<P, S) (2) 

with the covariant derivatives D^j^cj) = V f^cf) — ieiA^tj), D^^ = — ie2B^S, and the field strength tensors 
Ffii, — df_iA„ — d,yAf^, Hf^„ = d^Bi, — dvB^j, of the two U(l) gauge potential A^, i?^ with couphng constants 
ei and 62. The fields (f) and ^ are complex scalar fields (Higgs fields). 

We will study two different potentials, corresponding to Saffin 0] (Vi) and Rajantie et al. ^ (V2) 



vi{^,o = ^ (#* - vlf + Y (er - viy - A3 (#* - vl) (er - vl) (3) 

and 

V.i^^O = Y i'^^* ^llY + i^C - Vlf (4) 

In the first model^ (Vi) the interaction term makes the strings form a bound state depending on the 
value of A = A1A2 — 4A| for < A < ^y/XiM bounds states will be formed. In this paper we will 

investigate the special case where the parameters correspond to the Bogomol'nyi limit (more explicitly, 
to what would be the Bogomol'nyi limit if A3 = and G — Q) 

Al = 2el A2 = 2el (5) 

and we will eventually rescale parameters such that ?7i = ??2 = ei = 62 = 1, though we will keep them 
explicit throughout our equations. It is also worth reminding that by setting A3 = we recover two 
non-interacting Abelian Higgs strings in curved space-time. 

For the second case, the interaction ensures that it is energetically favourable to form bound states for 
Al > and A2 > (otherwise, the potential will not be bounded from below). This is clearly seen by 
observing that the potential energy of the ^ field vanishes where x = Oj therefore favouring that both 
string cores lie on top of each other. 



A. The Ansatz 

In the following we shall analyse the system of coupled differential equations associated with the 
gravitationally coupled system described above. This system will contain the Euler-Lagrange equations 
for the matter fields and the Einstein equations for the metric fields. In order to do that, let us write 
down the matter and gravitational fields as shown below. The most general, cylindrically symmetric line 
element invariant under boosts along the z— direction is: 

ds^ = N^{p)df - dp^ - L^{p)dip^ ~ N^{p)dz^ . (6) 



This type of potential has also been studied in relation to composite topological defects l2Cl l2l| . 
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The non- vanishing components of the Ricci tensor R'^ then read lj| : 



{LNN'y 2N" L" {N^L'Y 

where the prime denotes the derivative with respect to p. 



For the matter and gauge fields, we have 22l |: 



0(p,<^) = r?iMp)e™'^ , (8) 

^(P,^)=^]i/(p)e""'^, (9) 

A^da;'' = — (n - P{p))dip . (10) 
ei 

B.dx'" = — (m - R(p))dip . (11) 

62 

n and to are integers indexing the vorticity of the two Higgs fields around the z— axis. 

B. Equations of motion 

We define the following dimensionless variable and function: 

X = emip , L{x) = r]ieiL{p) . (12) 
Then, the total Lagrangian only depends on the following dimensionless coupling constants 

7 = SnGijl , g = — , g = — , = ^ , i = 1, 2, 3 or /3, = ^ , i = 1, 2 . (13) 

Varying the action with respect to the matter fields and metric functions, we obtain a system of six 
non- linear differential equations. The Euler-Lagrange equations for the matter field functions read: 

{N^_^ IdV^ .-12 (14) 
= 2g'fR, (17) 



V L 
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where the prime now and in the following denotes the derivative with respect to x. 
We use the Einstein equations in the following form: 



i?, 



(18) 



where T is the trace of the energy momentum tensor T = and the non-vanishing components of the 
energy- momentum tensor are (we use the notation of [l^ l with i — 1,2: 

(Tq )'i = es + ev + f w + Ui 
(T^)i = -es - ey + Eu, + Ui 

~ ~ 

(To). 



where 



and 



{h'r + iff 



[P'Y {R 



l\2 



2L2 2g2L2 



L2 L2 



-2 = |(/.^-ir+|/.M/^-.^r 

corresponding to the choice of potential Vi and V2 , respectively. 
We then obtain 



{LNN'Y 



and: 



= -7 



= 7 



(P') (ff 



l\2 



2L2 2g2i2 

2/t2p2 2i?2/2 (P')2 (i?')2 



L2 



2i2 ■ 2527:2 



, i = l,2 



(19) 



(20) 



(21) 



(22) 



(23) 



C. Boundary conditions 

The requirement of regularity at the origin leads to the following boundary conditions: 

h{0) = 0, /(O) = , P(0) = n , i?(0) = TO 

for the matter fields and 

iV(0) = 1, N'{0) = 0, L{0) = , i'(0) = 1 . 

for the metric fields. 

In the special case where n ^ and to = the boundary conditions change according to: 



(24) 



(25) 



h{Q) = 0, /'(O) = , P(0) = n , P(0) = . 



(26) 
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while for a n = 0, m 7^ string, they read 

h'{0) = 0, /(O) = , P(0) = , i?(0) = TO . (27) 

By letting the derivatives of the non-winding scalar field be zero at the origin instead of imposing the 
boundary conditions for the fields themselves, the non-winding scalar field can take non-zero values at 
the origin and, as it is indeed the case, a "condensate" of the non-winding scalar field will form in the 
core of the winding string. 

The finiteness of the energy per unit length requires: 

h{oo) = f , f{oo) = q , P(oo) = , R{oo) = . (28) 
We define as inertial mass per unit length of the solution 

E„, = J ^^T°dxdip (29) 
where gs is the determinant of the 2 + 1-dimensional space-time given by (t, (f). This then reads: 

poo 

E,n ^271 NL {es +e^ + + u,) (30) 



As mentioned earlier, for Baffin's model Vi, (3^ — corresponds to two non- interacting strings. In the 
BPS limit (the limit we are dealing with) the inertial mass (in units of 27r) of the (1, 1) string is just 
Ein{'^) = 2. This is because the mass is equal to the sum of the masses of the (l,0)-string and of the 
(0, l)-string, which is E^^^'^\j) — E^^^'^\'j) = 1. Correspondingly, the deficit angle (in units of 27r) is just 
6/(2'k) = 27. 

We can then define the binding energy of an (n, to) string as 

Ej^' = Et""^ ~ nEt'' - mE^'' (31) 



III. NUMERICAL RESULTS 
A. Generalities 

The solutions for 7 = 0, i.e. in flat space-time have been discussed in [5| and [6|. When 7 is increased 
from zero, the (p,q)-strings get deforme d by gravity. If 7 is not too large, globally regular solutions exist. 



These solutions fall into two categories [ij, ll5|: the string solutions and the Melvin solutions, where the 
latter have no flat-space counterpart. The string solutions are those of astrophysical interest since they 
describe solutions with a deflcit angle and because the Melvin solutions are heavier than their string 
counterparts. The metric functions then have the following behaviour at infinity: 

N{x 00) ci , L{x 00) C2X + C3 , C2 > , (32) 

where the Ci, i = 1, 2, 3 are constants that depend on the choice of 7, /3i, i = 1, 2, 3, <? and g. The solution 
has a deficit angle which is given by 



6 = 27r(I - C2) 



(33) 
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FIG. 1: The profiles of the metric functions A'^(a;) and L{x) (left) and of the matter functions P{x) = R(x), 
h{x) — f{x) (right) are shown for a Melvin solution (solid) and a string solution (dashed). Here Pi = f32 — 2, 
/Ja — 0.99, q — g — 1, J = 0.5 and n — m — 1. 

If the deficit angle S > 2n, C2 < and the solution would have a singularity at a finite, parameter- 
dependent value oi X = xq with L{x = xq) = 0, while N{xq) stays finite. These solutions are the so-called 
supermassive string solutions 13] (or inverted string solutions). 

The Melvin solutions exist for the same parameter values as the string solutions, but have a different 
asymptotic behaviour: 

N{x ^ cx)) -> aix^/^ , L{x oo) a-ix^'^l'^ , (34) 

where again oi and 02 are parameter dependent constants. This latter solution is thus not a cosmic string 
solution in the standard sense since it corresponds to a solution for which the circumference of circles will 
decay in the asymptotic region. 

The profiles of typical string and Melvin solutions are shown in Fig. [TJ 

B. String solutions for the potential Vi 

We have solved the above system of equations for /3i = 2, i — 1,2, g = 1 { which for 7 = corresponds 
to the BPS limit) and q = 1- First, we have fixed 7 = 0.5 and investigated the dependence of the deficit 
angle on the binding parameter /Sa and on the choice of n and m. Note that in order for h{x) — 1, 
f{x) = 1 to be the global minimum of the theory for our choice of the remaining parameters given above, 
we have to require < (3^ < 1. 

Om numerical results are given in Fig. [2l The deficit angle is decreasing for increasing (i^ and tends to 
zero in the limit /^s ^ 1. This is expected, since the increase of leads to a stronger binding between 
the strings, the total energy decreases and thus results in a smaller deficit angle. Moreover, the deficit 
angle for a fixed value of is smaller for an (n, n) string than for an (n + 1, n — 1) string, which itself 
is smaller than that for an {n + 2,n — 2) string etc. This is clearly seen for the deficit angle of the (2, 0) 
string in comparison to that of the (1, 1) string and accordingly for strings with n + m = 3 and n + m — 4, 
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FIG. 2: The value of the deficit angle S (in units of 2n) is shown as function of /Js for 7 — 0.5 and different values 
of (n, m). 

It is also obvious from Fig. [2] that globally regular string solutions exist only for a specific parameter 
range. For /Ja — 0, the deficit angle is given by (5/(27r) — {n + 'm)j. Since globally regular string solutions 
exist only for S < 2tt, we need {n + 771)7 ^ 1 or for 7 — 0.5, n + m < 2. Thus, for our specific choice 
of 7 and for (33 — we would not expect regular string solutions to exist for n + m > 3. This is clearly 
seen in the figure, where string solutions with deficit angle less than 27r exist for all values of only for 
the (1,0), (1, 1) and (2,0) cases. However, regular string solutions exist also for n + rn > 3 if P3 is large 
enough. For n + m > 3 and a fixed value of 7, regular string solutions exist for > I3'"^^\'^, n, m). We 
find e.g. /3^'"^^ (0.5, 2,1) w 0.76, /3^'=''^(0.5, 3, 0) « 0.86, /J^'^''^ (0.5, 2, 2) « 0.89 and /3^'=''^(0.5, 4, 0) « 0.94. 
The value (3^'^^'^ increases for n + m increasing, since a higher total winding requires a higher binding 
energy in order to render the string not too heavy. Moreover, for (n, m) — [n + k,n — k) strings, where 
< fc < n, the value of (3^'^^^ increases for increasing k. 

We have also studied the dependence of the deficit angle on the gravitational coupling for a fixed value 
of /33 — 0.95. The results are given in Fig. O 

As expected the deficit angle increases for increasing 7 and tends to 27r for 7 ^'-'^^\P3,n, m). Thus, 
globally regular strings exist only for 7 < j'^'^^^Ps, n, m). This critical value of the gravitational couphng 
decreases for n + m increasing and for a fixed value of n + m and (n, m) = {n + k,n — k), < k < n, 
decreases for k increasing. 

Clearly, the specific choice of /Ja, 7 and (n, m) determines whether globally regular solutions with deficit 
angle smaller than 27r exist. Increasing /Ja leads to a stronger binding between the strings, i.e. the energy 
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FIG. 3: The value of the deficit angle S (in units of 2n) is shown as function of 7 for fSg — 0.95 and different 
values of (n, m). 

of the (ri, m)-string is lowered and hence globally regular solutions exist for higher values of 7. This is 
demonstrated in Fig. |4] for (1,0), (1,1) and (2,1) strings. Here, we plot the value of the gravitational 
coupling up to which the regular solutions exist, 7^'^'"-', as function of fS^. Note that regular solutions exist 
only in the parameter domain below the respective curves. It is clear that the solutions exist for higher 
values of 7 when (3^ is increased. This is related to the increase in binding energy, e.g. for the (2, 1) 
solution with 7 = 0.5, we find that w -0.551 for f]^ = 0.8, while « -0.652 for /Jg = 0.95. In 

addition to the binding due to the potential, there is an additional binding due to gravity. The larger 7, 

(2 1) (2 1) 

the stronger is the binding: e.g. for = 0.9, we find E'^.J, w —0.621 for 7 — 0.5, while £'^j„ ~ —0.581 

for 7 = 0, in agreement with what it has been reported [5|. 

C. String solutions for the potential V2 

First, we have studied the solutions for j3i — (32 — ^, q ~ 9 ~ ^ and 7 = 0.5 to understand what 
qualitative differences appear in comparison to the potential Vi . While in the case of the Vi potential 
(and in the BPS limit) a (1, 0)-string is exactly equivalent to a (0, 1) string, this is different here. As in the 
Vi case, when one string is winding and the other is not, the non- winding string can form a condensate at 
the core of the other string. However, in this case, a condensate of the ^ field in the core of a </> string does 
not lower the energy; whereas a condensate of the </) inside the ^ core does. More specifically, for n = 1, 
m = 0, the scalar field function f{x) is not forced to zero at the origin and could develop a condensate. 
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FIG. 4: The value of i''^"' is shown as function of /Sa for (1, 0), (1, 1) and (2, 1) strings. Note that globally regular 
solutions for the specific choices of (n, m) exist only in the parameter domain below the corresponding curves. 



However, f{x) = q everywhere and correspondingly R{x) = gives the minimum for the potential for 
the boundary conditions we have. On the other hand, for n = and m = 1, the scalar field function h{x) 
is also not forced to zero at the origin, but does develop a condensate. At the origin, a value of h{x) < 1 
is the value that minimizes the energy best. Thus, a (l,0)-string is different from a (0, 1) string. This 
can also be seen when looking at the numerical values of the energies and deficit angles. For (1,0), we 
find E,r, « 0.859 and 5/{2tt) « 0.425, while for (0, 1), we have E,^ « 0.848 and S/{2tt) « 0.418 (see also 
Fig.©. 

In Fig. [5] we show the dependence of the deficit angle on the gravitational coupling 7. This is qualita- 
tively very similar to the case with the Vi potential. Again, we observe that globally regular string like 
solutions exist only up to a maximal value of the gravitational coupling 7 = 7^^''^(n, m). If one would 
increase 7 further, the solutions would have a deficit angle > 2t: and would thus be singular. 7''^''-'(n, to) 
decreases for n + m increasing and for a fixed value of n + to and (n, to) = (n + fc, n — fc), < fc < n, 
decreases for k increasing. 

As for the Vi case, the specific choice oi (3^, 7 and (n, m) determines, whether globally regular solutions 
exist. The fact that the strings interact and form bound states, lowers their energy. Therefore, for values 
of 7 and (n, m) where one would not expect regular global solution in the non-interacting case, regular 
solutions exists when strings interact. 
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FIG. 5: The value of the deficit angle S (in units of 2n) is shown as function of 7 for (3i = f32 = 1, q ~ g = 1 and 
different values of (n,m). 

D. Melvin solutions 

We have also studied the corresponding Melvin solutions for both potentials. Like their string coun- 
terparts, they exist only in a limited domain of the parameter plane. 

In the Vi case, for a fixed value of 7, P^'^^'' is equal for the string and Melvin solutions. This has already 
been observed in the case of a single gravitating string QQ. At this the inertial energies of 

the two types of solutions become equal. This is shown in Fig. [6] for (n, m) = (3, 0) and (n, m) = (2, 1), 
respectively. At the critical value of /Sg'^'^^ the curves of the inertial energies of the two types of solutions 
merge. Note also that for all cases studied in this paper, the inertial energy of the Melvin solutions is 
larger than that of the corresponding string solutions. 

For P3 < the string solutions become inverted/supermassive string solutions with a zero of 

the function L(x) at some parameter dependent value oi x = xo, while N{x = xq) stays finite. The 
Melvin solutions - on the other hand - become so-called Kasner solutions with N{x = xq) = and 
L(x — > xq) — !■ 00 at some finite and parameter-dependent value oi x — xq. The inverted string and 
Kasner solutions thus exist only on a finite range of the coordinate x and are thus sometimes also called 
"closed" solutions. These solutions have singularities in their metric, and it is out of the scope of this 
paper to investigate their implications. 

We have also studied the Melvin solutions in the case of V2. Again, the qualitative results are very 
similar to the case with Vi. In Fig. [71 we show the dependence of the inertial energy Ein on the 
gravitational coupling 7 for fii — 132 — i,q — g = ^ and three different choices of (n, m) . The inertial 
energy of the string solutions is always lower than that of the Melvin solutions and the two branches 
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FIG. 6: The dependence of the inertia! energy Ei„ on the binding parameter /Js is shown for string and Melvin 
solutions with (n, m) = (3, 0) (solid) and (n, m) = (2, 1) (dashed). Here Pi = P2 = 2, q = g = 1 and 7 = 0.5. 



of solutions meet at 7'-'^'^'', the maximal value of the gravitational coupling, beyond which no globally 
regular solutions exist. As before, the maximal value of the gravitational coupling decreases with n + m 
increasing. We include both (1,0) and (0,1) solutions, since, as stated above, the do differ for the V2 
case. 

It has been observed previously [l^ that the inertial energy of the Melvin solutions doesn't depend 
strongly on the winding number. The reason for this is that the main contribution to the energy of the 
Melvin solution comes from the gravitational field and is thus "insensitive" to the actual core structure 
of the solution. This explains why the energy of the Melvin solution practically doesn't depend on 
(see Fig. [5]), while it has a strong dependence on the gravitational coupling (see Fig. [7]). One could see 
these effects by studying the gravitationally active mass, the so-called Tolman mass Mtoi 23|, which for 
our definition of the metric is given by 14|: Mtoi oc limx^oc{LNN'). Clearly, for the standard string 
solutions, this mass is zero, while for the Melvin solutions it is non- vanishing. The Tolmann mass of the 
Melvin solutions would then depend on the windings of the strings (see also the discussions in ^15.]). 



IV. CONCLUSIONS AND DISCUSSION 

In this paper we have studied the gravitating properties of strings arising in field theoretical approaches 
to (p,q)-superstrings, by minimally coupling the field theoretical models to gravity. We find that there 
exist globally regular solutions of both the string type (solutions with conical deficit) and the Melvin 
type for a finite domain in parameter space. 

The existence of globally regular solutions depends very strongly on the (inertial) energy of the under- 
lying string specimen. Therefore, the specific choices of windings of the strings and the string interaction 



13 



2.6 




(1,0) 
(0,1) 

(1,1) 



E. 



in 



0.8 - 



0.6 







0.2 



0.4 



0.6 



0.8 



1.2 



1.4 



Y 



FIG. 7: The dependence of the inertia! energy Ei„ on the gravitational couphng 7 is shown for string and Melvin 
solutions in the case of V2 with I3i = — 1, q = g — 1. Here, we compare the cases {n,m) — (1,0) (solid), 
(n,m) = (0,1) (dashed) and {n,m) — (1,1) (dotted-dashed) . 

(as well as obviously the strength of the gravitational interaction) will determine whether regular solu- 
tions exists. Due to the interaction term, the inertial energy of a (p,q) string is lower than the sum of the 
inertial energy of its constituents; thus, the parameter region where regular strings exist gets enlarged 
with a stronger interaction term. We also found that it is not only the interaction between strings which 
makes the binding stronger, but a stronger gravitational interaction increases the binding further. 

Besides string-like solutions, we have also studied Melvin solution for these models. Like their string 
counterparts, they exist for a finite range in parameter space, and the range gets enlarged when the 
interaction between strings gets stronger. There is some critical value of the parameters beyond which 
no regular solution exists. At that critical value, the inertial energy of the Melvin and the string solution 
coincide, i.e. the solutions are degenerate. For regular solutions at values of the parameters which are 
not critical, the Melvin solution is always heavier than the string solution. 

While we have studied the situation where we have straight infinitely long strings, it would be very 
interesting to study the situation where three segments of strings meet in a Y junction. The gravitational 
lensing of such a Y junction has been studied for cosmic superstrings using several approximations [24]. 
It would be interesting to see whether we can learn something from the full study of these Y-junctions 
in field theory, and how they could be translated to proper cosmic superstrings. 

Another interesting situation would be the following: suppose that we have two strings with winding 
numbers such that their regular solution exists for a given set of parameters. Let us suppose also that 
the parameters are such that when these two strings interact to form a new bound string, there is no 
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regular solution of these new string. Then, we would have formed a Y junction in which one of the legs 
is a "supermassive string" with a zero of the metric at some finite value of the radial coordinate. This 
situation would be worth studying. 

There are other extensions of this work which deserve being studied. Even though we are dealing 
with "extensions" of the Abelian Higgs model here, string solutions in non-Abclian models have also 
been discussed 25] . It would be interesting to study the gravitational effects of the bound states in such 
models. 
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